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Abstract 

Based on general relativity, it can be argued that deviations from a uniform Hub- 
ble flow should be thought of as variations in the Universe's expansion velocity 
field, rather than being thought of as peculiar velocities with respect to a uniformly 
expanding space. The aim of this paper is to use the observed motions of galax- 
ies to map out variations in the Universe's expansion, and more importantly, to 
investigate whether real variations in the Hubble expansion are detectable given 
the observational uncertainties. All-sky maps of the observed variation in the ex- 
pansion are produced using measurements obtained along specific lines-of-sight and 
smearing them across the sky using a Gaussian profile. A map is produced for the 
final results of the HST Extragalactic Distance Scale Key Project for the Hubble 
constant, a comparison map is produced from a set of essentially independent data, 
and Monte Carlo techniques are used to analyse the statistical significance of the 
variation in the maps. A statistically significant difference in expansion rate of 9 
km s _1 Mpc -1 is found to occur across the sky. Comparing maps of the sky at dif- 
ferent distances appears to indicate two distinct sets of extrema with even stronger 
statistically significant variations. Within our supercluster, variations tend to occur 
near the supergalactic plane, and beyond our supercluster, variations tend to occur 
away from the supergalactic plane. Comparison with bulk flow studies shows some 
concordance, yet also suggests the bulk flow studies may suffer confusion, failing to 
discern the influence of multiple perturbations. 
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1 Introduction 



Conventionally, the Hubble flow is thought of as being completely uniform and 
isotropic. Deviations from a uniform Hubble flow are eliminated by imparting 
objects' observed residual recessional velocities into peculiar velocities, such 
that objects are thought to move with respect to a uniformly expanding space. 
However, empirically it is only valid to consider the velocity field of the matter 
and how everything is moving relative to everything else in the Universe. It is 
not possible to infer the existence of an absolute space that expands uniformly 
and that objects have peculiar velocities with respect to. Thus, deviations 
from a uniform Hubble flow should properly be considered deviations in the 
Universe's expansion itself. 



Interestingly, RaychaudhurJ (jl955 ) showed that (ignoring vorticity) if a ve- 



locity field has locally isotropic expansion, then the space is locally isotropic. 
Yet we know from examples such as gravitational lensing that inhomogeneities 
alter the curvature of space such that it is not locally isotropic. Thus, since 
space is not locally isotropic, then the Universe's expansion can not be locally 
isotropic either. Whether to conceive of the Universe expanding non-uniformly 
or whether to conceive of it expanding uniformly with superimposed peculiar 
velocities is more than just a conceptual issue, however. 



According to Raychaudhuri's equatiorj (jl955 ). the existence of shear in a ve 



locity field will lead to a decrease in the volume expansion. Since inhomo- 
geneities should introduce tidal forces and shear the velocity field, then the 
existence of overdensities and underdensities in the Universe should lead to 
shear throughout the Universe that decreases the Universe's volume expansion 
compared with that of a homogeneous universe. This effect should only be sig- 
nificant when measured locally in the vicinity of an inhomogeneity: the global 
influence should be quite small. Raychaudhuri's equation also shows that the 
existence of vorticity (and also velocity dispersion in the Newtonian version) 
will lead to an increase in the volume expansion. When structures start to 
collapse in the Universe and eventually become supported by vorticity or ve- 
locity dispersion, those regions of space cease shrinking, which can lead to 
an increase in the global expansion of the Universe. Thus, it is important to 
consider the influence inhomogeneities may have on the Universe's expansion. 

The Cosmological Principle — that the Universe is homogeneous and isotropic — 
is generally assumed to hold, since averaged over large enough scales the Uni- 
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verse will appear homogeneous. However, general relativity is needed to un- 
derstand not only small dense systems, but large diffuse systems such as the 
Universe, and according to Einstein's field equations, the spacetime corre- 
sponding to a homogeneous universe can not be used to represent a spatially- 
averaged inhomogeneous universe. This is because Einstein's field equations 
do not equate the spacetime to the mass-energy distribution directly. The 
energy-momentum tensor T a b depends on the Ricci tensor R a b and scalar R, 
which stem from taking derivatives of the metric tensor g^, with Einstein's 
equations equating 

R a b — -Rg a b = KT ab . 



If the left-hand side of the field equations for a homogeneous universe is 
equated to the spatially- averaged mass-energy of an inhomogeneous universe, 
there will generally be a discrepancy between the two sides of the field equa- 
tions, which will act like a cosmological constant and either accelerate or 
decelerate the universe's expansion from that expected for a homogeneous 
universe. Thus, even if the Universe may look homogeneous on large enough 
scales, assuming the Universe to expand uniformly is ultimately misleading. 
Several researchers have suggested this effect may even expla in the Universe's 
apparent acceleration (reported by Perlmutter et al. |, Il999 ) as b eing d ue to 



structure fo rmati on — Bildhauer and Futa masj ( 199ll ). Bene et alj (12003 ). and 



Kolb et alj (|2005j)— although Russ et alj (Il997f ) argue that the effect of inho- 



mogeneities should be small. 

Also, conceiving of the Universe's expansion as uniform and assigning the 
galaxi es peculiar velocities, bulk flow studies such as that of Hudson et alj 



(120041 ) have continued to find that the peculiar velocities with respect to the 
Cosmic Microwave Background (CMB) frame are correlated such that volumes 
of space of order 100 Mpc in radius are moving with bulk velocities of approx- 
imately 300-700 km s _1 . This suggests inhomogeneities significantly perturb 
the velocity field of the Universe. The existence of the Universe's large-scale 
structure of voids and superclusters suggests the voids are underdense regions 
that have been decelerated less due to gravity so they have ballooned up into 
roughly spherical regions without undergoing structure formation, while the 
superclusters are overdense regions where gravity has overcome the Universe's 
expansion such that they have reached turnaround and collapsed in their dens- 
est regions. 



Moffat and Tatarskj ( 1995 ) looked at what observational effects we would theo- 



retically observe if we were to inhabit a local void. Via comparison of their the- 
oretical curves with a survey of redshift-distance determinations, they found 
the data were better fit by a model with a local void than by a homogeneous 
universe. Zehavi et alj ( 1998 ) used 44 type la supernova H values to show 
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that we may just inhabit an underdense region of the Universe (where the ex- 
pansion in the velocity field has been slowed less due to gravity than i n mor e 
dense regions of the Universe). Referring to fig. 4 of Freedman et alj ( 2001 ). 
it appears that the H values tend to fall off beyond a distance of 100 Mpc, 
which suggests the Universe may be expanding faster locally. A here-there 
difference in the Universe's expansion could be an alternative to the notion of 
a now-then difference, which is th e assumption the Universe's supposed accel- 
eration (Perlmutter et al. |. Il999 ) rests on, so it is important to account for 
the possible influence of inhomogeneities on the Universe's expansion if the 
cosmological parameters are to be properly determined. 



Thus, in this paper we will not assume the existence of a uniform spatial 
expansion with peculiar velocities superimposed. We will use Hq values mea- 
sured along different lines-of-sight to see whether local variation in Ho exists, 
and to produce all-sky maps of the observed variation across the sky. If more 
variation exists in the maps than should be expected due to measurement er- 
rors in the data, and if the high and low values of H are correlated in position 
on the sky, then this will be taken as evidence that the expansion is indeed 
locally anisotropic across the sky. Since bulk flow studies find bulk flows of a 
few hundred km s _1 on 100 Mpc scal es, wh ich is predicted depending on the 
cosmological model (e.g. see Zaroubi L l2002h . and bulk flows only show the net 
flow of a sample volume rather than the individual variations in the velocity 
field, then it would be expected that variations in H observed on this scale 
should be at least a few km s -1 Mpc -1 . 



While it is easy enough to measure how fast objects are expanding away 
from us via redshifts, it is the determination of accurate distances that is 
problematic in the determination of Ho- Historically, the errors in Hq have 
been so great that it would be difficult to study real variation in the Universe's 
expansion rate. The most accurate work to date to stu dy H p is the HST 
Extragalactic Distance Scale Key Project (Freedman et al. 1. 12001 . hereafter the 
HST Key Project), which yielded distances accurate enough for a meaningful 
study of real variation in Ho, especially since most of the errors are systematic 
and shared by all the H determinations. Thus, we will map the directional 
variation in Hq using the HST Key Project data, comparing with a second 
set of data to examine whether the same general trend is observed. The data 
selection will be discussed in Sect. 2. The technique used to generate all-sky Hq 
maps and study the significance of the variations and the impact of distance 
will be outlined in Sect. 3. In Sect. 4 the results will be examined from various 
frames of reference and a comparison with bulk flow studies will be made. 
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2 Selection of Data 



2.1 HST Key Project data 



The HST Key Project data for H have been selected as the primary data set. 
This set offers a reasonably large set of values that is distributed nicely for 
the purpose of making all-sky maps. Also, despite the fact that the H values 
depend on the data of other researchers, the values were all analysed by the 
HST Key Project team to be consistent with each other. Thus, the systematic 
errors should be similar in most cases to minimize the effect on the study of 
variations in H so that most of the uncertainty in the relative values of H 
will just be in the random errors. 



There are 74 value s that were published in the final HST Key Project paper 
( Freedm an et al. I, l200lh as well as 2 values that were obtained from Sakai 
| ( 2001 ). The data stem from the use of 4 different methods to obtain H : 
the Tully-Fisher relation, the fundamental plane relation, surface brightness 
fluctuations, and type la supernovae. 



In all cases, the H values used are those for which the recessional velocities 
have been corrected to the CMB frame. That way the Hq values are all in one 
frame of reference, rather than making corrections for infalls, which would 
be contrary to our purposes. Thus, the most recent distances from Freedman 
et al. and Sakai were used, and where H values corresponding to non-CMB 
recessional velocities were reported, t he CM B-frame recession al velo cities from 
the former papers of Ferrarese et alj fl2000h and Sakai et alj fl2000f ) were used 
to yield the CMB-frame Hq values. 



The 2 value s obtained from Sakai (j2001 ) had originally appeared in Sakai 



etalJ (l2000h but did not appear with modified distances in Freedman et alj 



One of these values was for Pavo 2, which is a component that was 



separated from the Pavo Cluster due to some of the galaxies yielding different 
Tully-Fisher distances, but since the galaxies had similar recessional velocities, 
they yielded quite different H values. For our purposes, it only seems fair to 
include the values for both Pavo components. 



Reported in Table 1 are the objects used, their celestial co-ordinates, and their 
Hq values with corresponding 1—a random errors. The celestial co-ordinates 
for the objects were obtained via SIMBAD. 



5 



Table 1 

HST Key Project H Q 



Data 



Object ID 
or Cluster 



a 

(hours) 



S 

(degrees) 



Ho 

(km/s/Mpc) 



Reference 



SN 1990O 
SN 1990T 
SN 1990af 
SN 1991S 
SN 1991U 
SN 1991ag 
SN 1992J 
SN 1992P 
SN 1992ae 
SN 1992ag 
SN 1992al 
SN 1992aq 
SN 1992au 
SN 1992bc 
SN 1992bg 
SN 1992bh 
SN 1992bk 
SN 1992bl 
SN 1992bo 
SN 1992bp 
SN 1992br 
SN 1992bs 
SN 1993B 
SN 19930 
SN 1993ag 
SN 1993ah 
SN 1993ac 
SN 1993ae 
SN 1994M 



17.15 
19.59 
21.35 
10.29 
13.23 
20.00 
10.09 
12.42 
21.28 
13.24 
20.46 
23.04 
00.10 
03.05 
07.42 
04.59 
03.43 
23.15 
01.22 
03.36 
01.45 
03.29 
10.35 
13.31 
10.03 
23.52 
05.46 
01.29 
12.31 



+16.2 
-56.2 
-62.4 
+22.0 
-26.1 
-55.2 
-26.4 
+10.2 
-61.3 
-23.5 
-51.2 
-37.2 
-49.6 
-39.3 
-62.3 
-58.5 
-53.4 
-44.4 
-34.1 
-18.2 
-56.1 
-37.2 
-34.3 
-33.1 
-35.3 
-27.6 
+63.2 
-01.6 
+00.4 



67.3 ± 2.3 

75.6 ± 3.1 
75.8 ± 2.8 
69.8 ± 2.8 

83.7 ± 3.4 

73.7 ± 2.9 

74.5 ± 3.1 

64.8 ± 2.2 

81.6 ± 3.4 

76.1 ± 2.7 
72.8 ± 2.4 

64.7 ± 2.4 

69.4 ± 2.9 
67.0 ± 2.1 

70.6 ± 2.4 

66.7 ± 2.3 

73.6 ± 2.6 

72.7 ± 2.6 

69.7 ± 2.4 
76.3 ± 2.6 

67.2 ± 3.1 

67.8 ± 2.8 

69.8 ± 2.4 

65.9 ± 2.1 
69.6 ± 2.4 
71.9 ± 2.9 
72.9 ± 2.7 
75.6 ± 3.1 
74.9 ± 2.6 



Freed. . f2001 
Freed. 
Freed. 
Freed. 
Freed. 
Freed. 
Freed. 
Freed. 
Freed. 
Freed. 
Freed. 
Freed. . (2001 
Freed. . C2001 
Freed. 
Freed. 
Freed. 
Freed. 
Freed. 
Freed. 
Freed. 
Freed. 
Freed. 
Freed. 
Freed. 
Freed. 
Freed. 
Freed. 
Freed. 
Freed. 
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Table 1 
- continued 
Object ID 

or Cluster 



a 

(hours) 



8 

(degrees) 



H 

(km/s/Mpc) 



Reference 



SN 1994Q 


16 


50 


+40 


3 


SN 1994S 


12 


31 


+29 


1 


SN 1994T 


13 


21 


-02 


1 


SN 1995ac 


22 


45 


-08 


5 


SN 1995ak 


02 


45 


+03 


1 


SN 1996C 


13 


50 


+49 


2 


SN 1996bl 


00 


36 


+11 


2 


Abell 1367 


11 


74 


+19 


8 


Abell 2197 


16 


47 


+40 


9 


Abell 262 


01 


88 


+36 


1 


Abell 2634 


23 


64 


+27 





Abell 3574 


13 


82 


-30 


3 


Abell 400 


02 


96 


+06 


6 


Antlia 


10 


50 


-35 


3 


Cancer 


08 


35 


+21 





Cen 30 


12 


77 


-41 





Cen 45 


12 


80 


-40 


4 


Coma 


13 


00 


+28 





Eridanus 


00 


50 


-21 


5 


ESO 508 


13 


17 


-23 


1 


Fornax 


03 


64 


-35 


5 


Hydra 


10 


61 


-27 


5 


MDL 59 


22 


01 


-32 


2 


NGC 3557 


11 


17 


-37 


5 


NGC 383 


01 


12 


+32 


4 


NGC 507 


01 


39 


+33 


3 


Pavo 


21 


23 


-57 


8 


Pavo 2 


21 


23 


-57 


8 


Pegasus 


23 


34 


+08 


2 



68.0 ± 2.7 
72.5 ± 2.5 
71.5 ± 2.6 

78.8 ± 2.7 

80.9 ± 2.8 
66.3 ± 2.5 
78.7 ± 2.7 
75.2 ± 12.5 
77.2 ± 12.5 
70.9 ± 11.8 

77.7 ± 12.4 

76.2 ± 12.2 

79.3 ± 12.6 

68.8 ± 11.3 

67.1 ± 11.0 
75.8 ± 12.8 

70.7 ± 11.9 

83.5 ± 13.4 

77.6 ± 12.9 

79.8 ± 13.0 

92.2 ± 15.3 
69.6 ± 11.1 
73.6 ± 11.8 
85.0 ± 14.4 

73.9 ± 11.9 
84.9 ± 13.5 
124.4 ± 19.0 

86.3 ± 14.2 

66.4 ± 10.7 



Freed. _ (2001 
Freed. 
Freed. 
Freed. 
Freed. 
Freed. 
Freed. 
Freed. 

Sak. _ 
Freed. 
Freed. 
Freed. 
Freed. 
Freed. 
Freed. 
Freed. 
Freed. 
Freed. 
Freed. 
Freed. 
Freed. 
Freed. 
Freed. 
Freed. 
Freed. 
Freed. 
Freed. 

Sak. . (2001) 
Freed. . (2001 ) 
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Table 1 
- continued 
Object ID 

or Cluster 



Q: 

(hours) 



(degrees) 



(km/s/Mpc) 



Reference 



Ursa Major 
Dorado 
Hydra I 
Abell S753 
GRM 15 
Abell 3574 
Abell 194 
Abell S639 

Coma 
DC 2345-28 
Abell 539 
Abell 3381 
NGC 4373 
Abell 262 
Abell 3560 
Abell 3565 
Abell 3742 
Coma 



11.95 
04.27 
10.61 
14.06 
20.05 
13.82 
01.53 
10.68 
13.00 
23.75 
05.28 
06.17 
12.42 
01.88 
13.53 
13.61 
21.11 
13.00 



+49.3 
-55.6 
-27.5 
-34.0 
-55.5 
-30.3 
-01.5 
-46.2 
+28.0 
-28.0 
+06.5 
-33.6 
-39.8 
+36.1 
-33.2 
-34.0 
-47.1 
+28.0 



54.8 ± 8.6 

81.9 ± 8.7 

82.8 ± 8.4 

87.5 ± 7.9 

95.6 ± 10.0 
92.0 ± 10.0 
91.3 ± 7.5 
109.7 ± 9.9 
83.2 ± 6.0 
83.2 ± 6.4 
86.2 ± 6.5 

88.9 ± 8.3 
99.9 ± 11.2 

69.0 ± 7.7 

78.1 ± 8.7 

70.2 ± 7.8 
70.0 ± 7.8 

70.3 ± 17.9 



Freed. (2001 



Freed. 
Freed. 
Freed. 
Freed. 
Freed. 
Freed. 
Freed. 
Freed. 
Freed. 
Freed. 
Freed. 
Freed. 
Freed. 
Freed. 
Freed. 
Freed. 
Freed. 



(2001 
(2001 
(2001 
(2001 
(2001 
(2001 
(2001 
(2001 
(2001 
(2001 
(2001 
(2001 
(2001 
(2001 
(2001 
(2001 
(2001 



2.2 Comparison data 



In an effort to test whether variations detected in the HST Key Project data 
exist independently of this data set, another set of data has been compiled. 
This data set has been constructed by conducting a literature search for de- 
terminations of Hq and using papers that report measurements of Hq along 
individual lines-of-sight, include uncertainties, and for which the observed ob- 
jects are within the distance range of the HST Key Project data. This yields 
57 values, which will hereafter be referred to as the comparison values. Re- 
ported in Table 2 are the objects used, celestial co-ordinates (obtained via 
SIMBAD), and Hq values with their reported errors. 

It should be noted that the comparison values stem from a mixed bag of 
methods and different researcher analysis; thus, the potential for this set to 
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Table 2 

Comparison Hq Data 



Object ID 
or Cluster 



a 

(hours) 



(degrees) 



(km/s/Mpc) 



Reference 



Coma 

Coma 

Virgo 
NGC 7331 

Virgo 

Fornax 
NGC 5846 
NGC 1365 

Coma 
IC 4051 

Coma 
NGC 4889 
NGC 3309 
NGC 4881 
Abell 2256 

Coma 
Abell 262 
Abell 3560 
Abell 3565 
Abell 3742 

Coma 
PGC 14638 
PGC 39724 
PGC 51233 
PGC 00218 
PGC 10208 
PGC 35164 
PGC 43798 

Fornax 



13.00 
13.00 
12.50 
22.62 
12.50 
03.64 
15.11 
03.56 
13.00 
13.00 
13.00 
13.00 
10.61 
13.00 
17.06 
13.00 
01.88 
13.53 
13.61 
21.11 
13.00 
04.20 
12.33 
14.34 
00.05 
02.70 
11.44 
12.89 
03.64 



+28.0 
+28.0 
+13.2 
+34.4 
+13.2 
-35.5 
+01.6 
-36.1 
+28.0 
+28.0 
+28.0 
+28.0 
-27.5 
+28.0 
+78.7 
+28.0 
+36.1 
-33.2 
-34.0 
-47.1 
+28.0 
-32.9 
+29.6 
+03.9 
+16.1 
+00.4 
+43.6 
+02.2 
-35.5 



71 ± 30 
78 ± 11 
80 ± 16 
70 ± 14 
87 + 7 

73 + 5 
65 + 8 
75 + 5 
75 + 6 
68 + 6 

70 + 7 

85 ± 10 
46 + 5 

71 ± 11 

72 ± 22 
67 ± 26 

82 + 8 

86 + 7 

83 + 6 
78 + 6 
60 ± 11 
60 ± 15 
44 ± 15 
60 ± 15 

58 ± 15 
60 ± 15 

59 ± 15 
39 ± 15 

74 + 5 



Her. . (1995 ) 
Whit. . (1995) 
Zas. . (1996) 



Zas. 

Ford 

Ford 

Forb. 

Mad. 

Gregg 

Baum 



1996) 

1996 ^ 

1996 ) 

(1996) 

(1996) 

(1997) 

(1997) 



Hjor. . (1997) 
Jen. . (1997) 
Jen. . (1997) 
Thorn. . (1997) 



Myers 
Myers 
Lauer 
Lauer 
Lauer 
Lauer 
Sal. (1998) 



1997) 
1997) 
1998) 
1998 ) 
1998 ) 



Pat. 
Pat. 
Pat. 
Pat. 
Pat. 
Pat. 
Pat. 
Tul. 



(1998 ) 
(1998) 
(1998 ) 
(1998 ) 
(1998) 
(1998 ) 
(1998 ) 
(2000) 
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Table 2 
- continued 
Object ID 

or Cluster 



a 

(hours) 



5 

(degrees) 



(km/s/Mpc) 



Reference 



Ursa Major 
Pisces Fil. 

Coma 
Abell 1367 
Antlia 
Cen 30 
Pegasus 
Hydra I 
Cancer 
Abell 400 
Abell 2634 
Abell 262 
Abell 496 
Abell 779 
Abell 1060 
Abell 1656(a) 
Abell 1656(b) 
Abell 2199 
Abell 2666 
Abell 3389 
Abell 3565 
Abell 3581 
Abell 3656 
Abell 3742 
NGC 4073 
NGC 4709 
NGC 5193 
Coma 



11.50 
01.12 
13.00 
11.74 
10.50 
12.77 
23.34 
10.61 
08.35 
02.96 
23.64 
01.88 
04.56 
09.33 
10.61 
12.99 
12.99 
16.48 
23.85 
06.36 
13.61 
14.12 
19.98 
21.11 
12.07 
12.83 
13.53 
13.00 



+55.0 
+32.4 
+28.0 
+19.8 
-35.3 
-41.0 
+08.2 
-27.5 
+21.0 
+06.6 
+27.0 
+36.1 
-13.2 
+33.8 
-27.5 
+28.0 
+28.0 
+39.6 
+27.1 
-65.0 
-34.0 
-27.2 
-38.3 
-47.1 
+01.9 
-41.4 
-33.2 
+28.0 



59 + 6 

79 + 2 
83 + 2 
77 + 2 
86 + 3 
83 + 3 
77 + 3 
70 + 2 

80 + 2 

76 + 2 

70 + 2 

77 + 4 
74 + 2 

73 + 2 

74 + 4 
79 + 3 
82 + 3 

71 + 2 
68 + 2 

70 + 2 

78 + 4 
74 + 2 

72 + 3 

81 + 4 
64 + 2 
107 ± 7 
85 + 6 

71 + 8 



Tul. . (2000) 
Tul. 
Tul. 
Tul. 
Tul. 
Tul. 
Tul. 
Tul. 
Tul. 
Tul. 
Tul. 
Jen. , 
Jen. _ 
Jen. . 
Jen. , 
Jen. _ 
Jen. . 
Jen. , 
Jen. _ 
Jen. . 
Jen. , 
Jen. _ 
Jen. i 
Jen. , 
Jen. _ 
Jen. i 
Jen. . 
Liu 
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be corrupted by systematics between values is greater than in the HST Key 
Project set. Also, various infall corrections have been made, and these are not 
always clearly stated in the papers, so these values may not consistently be 
in the CMB frame. One of the comparison values actually stems from a HST 
Key Project Cepheid distance value, while 4 of the HST Key Project values 
depend on comparison data surface brightness fluctuation distances, but the 
data sets are essentially independent. 



3 Contour mapping 

3. 1 Technique 

A method is required to generate an all-sky map based on a set of values 
located at specific positions on the sky. One method would be to fit spherical 
harmonics; however, this would require more higher order terms than could 
ever be convenient in order to not force structure into the map from the lower 
order terms. Our chosen method is to smear the H values over the sky using 
a Gaussian profile for each data point. 

The Gaussian smearing method involves laying out a grid on the sky and 
calculating weighted mean values of H at each grid point, weighting each 
actual data point in the average according to its angular separation 9 from the 
grid point such that the weightings fall off as a Gaussian. Thus, the weighting 
W of each data point is given by 

W = -^e- e2 /^ 

with the standard deviation a controlling how broad the smearing is. Contours 
of constant H are then interpolated within the grid of averaged Hq values to 
generate contour maps of H . The values are weighted only by their separa- 
tions, not their uncertainties (the type la supernova values are more distant 
and have smaller uncertainties, but the effect of distance on the map will be 
specifically explored in Sect. 3.3). 

Gaussian smearing succeeds in creating an all-sky map from a sample of data 
points associated with specific positions on the sky, and it also averages out 
the impact of errors associated with individual data points so that variations 
correlated with directions on the sky have the opportunity to manifest them- 
selves. This will be sufficient for studying large-scale variations in H , although 
there is no hope of studying any variations that do not have a large angular 
extent, as there is insufficient sky coverage. It should be kept in mind that 
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this method will also smear out the extrema for any actual variation though, 
so the range of any real variation in the maps will be diminished somewhat. 

For the contour maps presented in this paper, while the angular separations 
are calculated in spherical co-ordinates (using the dot product of unit vectors 
for the angular positions of the grid points and data points), the grid points 
are positioned for a cylindrical projection of the sky and are then mapped out 
in the form of a sinusoidal projection. The sinusoidal projection is a pseudo- 
cylindrical projection that preserves areas by keeping latitude lines parallel 
but shortening their length longitudinally according to the sine of the polar 
angle (or cosine of the declination). 

Unless stated otherwise, the grid points are set 1° apart in right ascension and 
1° apart in declination, as this appears to be a fine enough grid for the purpose 
of interpolating contours, and the Gaussian weighting profile falls off with a 
standard deviation of 25°, as this is approximately the typical separation of 
the real data points and is a sufficiently broad smearing to fill in the holes in 
the distribution of data on the sky. 

In Fig. 1 contour maps appear (in Galactic co-ordinates) for the HST Key 
Project data and the comparison data. The extrema are in similar directions 
for the two data sets and are similar in magnitude, except that one of the 
maxima in the comparison map is weaker. While the comparison map may 
not serve as a completely definitive cross-check of the reality of the variations 
observed in the HST Key Project map, it certainly seems to show some agree- 
ment, and it does not differ from the HST Key Project map as much as would 
be expected if the variation in Hq were due to uncertainties in the grid Hq 
values alone. This suggests the variation in the maps exists independently of 
the uncertainties in the particular determinations of H . Since there is more 
uncertainty in the comparison map than the Key Project map, the differences 
between the maps probably mostly reflect errors in the comparison map rather 
than the Key Project map. 

If the Pavo, Pavo 2, and Ursa Major determinations are removed from the 
HST Key Project data (due to their discrepant values), the map looks largely 
unchanged from the original: the primary maximum goes down slightly to 80 
km s -1 Mpc -1 and the primary minimum goes up slightly to 71 km s" 1 Mpc -1 
while the secondary extrema remain the same, so that the secondary extrema 
now become the primary extrema in this case. If PGC 39724 and NGC 4709 
are removed from the comparison data, again the map is largely unchanged: 
the primary maximum is slightly lower at just under 79 km s _1 Mpc -1 , while 
the other extrema appear unchanged. Looking at Fig. 1, it is apparent that 
the extrema tend not to be centred on the lines-of-sight to the individual data 
determinations, so the extrema are resulting from trends in the data, rather 
than from specific high or low H values. 
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Fig. 1. Hubble constant contour maps (in Galactic co-ordinates for a sinusoidal pro- 
jection of the sky) for (a) the 76 HST Key Project Hq values and (b) 57 comparison 
Ho values. Positions of the actual data points are indicated by triangles, and the 
contours range from low (dark) to high (light) values of Ho (in km s _1 Mpc^ 1 ) as 
indicated. 

3.2 Magnitude of variation 

As was previously discussed in Sect. 3.1, the Gaussian smearing lessens the 
range of any real variation. Using Gaussians with successively smaller standard 
deviations of 20°, 15°, and 10°, as the standard deviation gets smaller, the 
extrema are picked out with less smearing, but errors in the data also start 
to have a greater impact. At 10°, the range of variation is greater than 30 km 
s _1 Mpc -1 , but the grid values are mostly being determined by individual H 
values so the errors in the grid values approach those of the H determinations. 
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The actual extrema are separated by 38 km s _1 Mpc -1 , which may be an 
underestimate to the magnitude of any real variation if the standard deviation 
could approach zero, but the errors have already become so large that they 
are likely making the range of variation appear artificially large as it is. Being 
conservative, the range of variation appears to be ~30 km s" 1 Mpc -1 . Thus, 
large and small values for the standard deviation yield complementary aspects 
of the map: smaller standard deviations smooth out the range of variation less, 
while larger standard deviations yield the overall trend in the map with less 
error. 

One may question whether the variation is statistically significant or whether 
it is likely that this much variation would be found in the map due to measure- 
ment errors in H alone. Assuming that the 1-a random errors are accurate, 
the significance of the result can be tested statistically by using Monte Carlo 
simulations of the data. This is accomplished using two different methods: 
one which assumes the systematic differences between the mean H values 
are real for each of the 4 methods used to derive H Q , and one which assumes 
that the systematic differences stem from variations in H or errors. The first 
method involves calculating a new value of H Q for each position by using the 
weighted mean value of H corresponding to the actual distance method that 
was used to obtain the real H value, and then adding Gaussian deviates (var- 
ied over a 3-cr range) to the data using the weighted mean 1-a error for the 
corresponding method to simulate the expected scatter. The second method 
involves using the weighted mean value of H for all the data for the H value 
at each position and adding Gaussian deviates to the values according to the 
actual 1-a errors for each data point. 

For 10,000 different sets of simulated data for each Monte Carlo method, H 
maps are calculated and the differences between the maximum and minimum 
H Q grid point values are found. Comparing the random sky variations with 
the observed variations depends strongly on the value of the standard devia- 
tion used to smear the values over the sky. This is because smaller standard 
deviations smooth out the errors less and make it more likely to find greater 
variation in the maps. For the first method, the results are that the magnitude 
of variation is only as great as in the real data 5.61% of the time for a standard 
deviation of 45° (which yields a range of 6.3 km s" 1 Mpc" 1 ), 9.90% of the time 
for a standard deviation of 35° (which yields a range of 8.8 km s" 1 Mpc -1 ), 
or 37.47% of the time for a standard deviation of 25° (which yields a range 
of 12.9 km s _1 Mpc" 1 ). The results for the second method are respectively 
3.82%, 4.93%, or 13.08%. Thus, assuming the realistic case is somewhere be- 
tween the two methods, and assuming anything with a less than 5% chance is 
"statistically significant," it appears that a statistically significant difference 
of 6 to 9 km s -1 Mpc -1 can be demonstrated with the broader values of the 
standard deviation. 
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While the above Monte Carlo methods should give a reasonable measure of the 
statistical significance of the variation, they depend on the \-a uncertainties 
reported for the data being accurate. Independent of the error bars, one can 
study whether the values are correlated in position on the sky, with higher 
values tending to be near higher values and lower values tending to be near 
lower values. Randomly reassigning the H values to the actual H positions 
on the sky and computing maps for several randomizations reveals whether 
the actual map has more variation than just the scatter in the data should 
produce. If there is real variation in the data, this method will not yield accu- 
rate measures of the statistical significance, since there will be extra scatter in 
the data that will tend to allow more variation in the map than errors alone 
should produce. However, this method at least yields an upper limit to the 
likelihood that as much variation could occur in the map if there were no real 
variation in the data. 

For 10,000 randomizations of the above method, the results are that the range 
of variation is only as great as for the real map 14.78% of the time for a 
standard deviation of 45°, 21.55% of the time for a standard deviation of 35°, 
or 42.10% of the time for a standard deviation of 25°. Thus, even if one is a 
complete skeptic that there is no real variation in the data, there appears to 
be more variation than there should be, although of weaker significance than 
the above Monte Carlo methods found. On the other hand, if there is real 
variation in the data, these are upper bounds on the percentages, so this weaker 
significance is not inconsistent with the results of the above methods. Since 
this method does not depend on uncertainties, it can also be applied to the 
comparison data. For the complete set of 133 data points, the results are that 
the variation is only as great 1.16% of the time for a standard deviation of 45° 
(which yields a range of 6.6 km s -1 Mpc -1 ), 2.60% of the time for a standard 
deviation of 35° (which yields a range of 8.9 km s -1 Mpc -1 ), or 19.29% of 
the time for a standard deviation of 25° (which yields a range of 12.6 km s -1 
Mpc -1 ). Thus, the difficulty in producing coincidental correlations in larger 
sets of data allows the conglomerate data set to substantiate a statistically 
significant variation of ~9 km s -1 Mpc -1 . 

3.3 Directional uncertainty and distance dependence 

The directional uncertainty of the extrema in the map can be tested by adding 
Gaussian deviates (in a 3-<r range) to the 76 data points, computing the H 
grid map, and seeing how much this affects the positions of the maximum 
and minimum Ho grid values in the map. The extrema for each of 500 ran- 
domizations are plotted in Fig. 2. Each randomization is calculated with grid 
separations of 0.5° so that the extrema can be plotted with less granularity. It 
can be seen that the probability distribution for the positions of the extrema 
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Fig. 2. Hubble constant contour map of Fig. 1(a) with 500 random extrema for maps 
calculated with Gaussian deviates for the 76 HST Key Project Hq values. Positions 
of the minima and maxima are indicated by dark and light dots respectively, and 
the contours range from low (dark) to high (light) values of Hq (in km s -1 Mpc -1 ) 
as indicated. 

is not the same in all directions, but primary and secondary extrema exist 
with directional uncertainties of order 10° to 20°. 

To test the influence of depth on the map, an additional weighting factor for 
distance is added in the computation of the grid maps. Maps are computed 
for various nominal distances by weighting the data according to the fraction 
of the distance that is shared in common with a given nominal distance. For 
distances more distant than the nominal distance, the weightings go as the 
ratio of nominal distance to object distance. For distances closer than the 
nominal distance, the weightings go as the ratio of object distance to nominal 
distance. 

About half of the HST Key Project data are for objects between 13 and 70 
Mpc, with the other half between 70 and 467 Mpc. In Fig. 3 distance-weighted 
maps appear for 6 nominal distances: 30 Mpc, 50 Mpc, 80 Mpc, 120 Mpc, 180 
Mpc, and 300 Mpc. It is apparent that nearby, one pair of extrema from 
Fig. 1 dominates, and the magnitude of variation of this pair falls off with 
distance. The minimum is near (cv = 9 h 30 m , 5 = +70°), and the maximum 
is near (a = 19 ft 30 m , 5 = —70°). Less apparent is that the grid values of H 
at the secondary extrema from Fig. 1 remain roughly constant with distance 
and dominate only at the greatest distances where the range of variation from 
the first pair of extrema has become small enough. The secondary minimum is 
near (cv = 18 h m , 5 = +15°), and the secondary maximum is near (a = 5 h 30 m , 
5 = +5°). 
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Fig. 3. Hubble constant contour maps (in Galactic co-ordinates for a sinusoidal 
projection of the sky) for the 76 HST Key Project Hq values weighted for nominal 
distances of (a) 30 Mpc, (b) 50 Mpc, (c) 80 Mpc, (d) 120 Mpc, (e) 180 Mpc, and 
(f) 300 Mpc. Positions of the actual data points are indicated by triangles, and the 
contours range from low (dark) to high (light) values of Hq (in km s _1 Mpc -1 ) as 
indicated. 

It should be noted that more distant values of H will always sample local 
expansion as well, and while the converse is not true, the weighting method 
allows values of Hq for the full range of distances to affect all the maps to 
some degree. The extrema that dominate locally appear to imply some sort of 
local effect, which could be interpreted as a dipole due to a bulk flow of a local 
sample volume with respect to the CMB frame, since the extrema are roughly 
opposite on the sky. The secondary extrema remain constant with distance, 
and are also roughly opposite of each other on the sky, perhaps suggesting an 
independent bulk flow of a larger-scale volume. 

Removing all the data points from the HST Key Project data set that have 
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an uncertainty greater than 8 km s" 1 Mpc" 1 and producing a map with the 
remaining 44 points yields a map quite similar to the most distant maps of Fig. 
3. The remaining data essentially consist of type la supernova measurements 
and a few other surface brightness fluctuation and fundamental plane mea- 
surements that are for distances about 50 Mpc and greater. The map ranges 
from 68 km s" 1 Mpc" 1 at the minimum to 81 s" 1 Mpc -1 at the maximum, 
and the primary extrema of Fig. 1 now appear even weaker than in Fig. 3(f). 
Since these data mostly have uncertainties of less than 3 km s" 1 Mpc" 1 , this 
provides stronger support for the reality of the extrema observed on large 
scales in Fig. 3. 

It is difficult to make an all-sky map from only 76 points: weighting for dis- 
tance only makes it more difficult to make meaningful all-sky maps. Thus, 
to test the significance of the observed variation with distance, the HST Key 
Project data are combined with the comparison data, and then divided into 
two sets according to distance. Randomly reassigning the Hq values to the 
data positions for the nearby group of 67 data points (to get a limit on the 
statistical significance) yields the result that as much variation occurs only 
0.16% of the time for a standard deviation of 35° (which yields a range of 
variation of 18.8 km s" 1 Mpc" 1 ) or 7.20% of the time for a standard deviation 
of 25° (which yields a range of variation of 25.3 km s" 1 Mpc" 1 ). Likewise, the 
more distant group of 66 data points gives respectively 24.39% (corresponding 
to a range of variation of 6.0 km s" 1 Mpc" 1 ) or 23.05% (corresponding to a 
range of 11.1 km s" 1 Mpc" 1 ). Thus, given that these percentages are upper 
limits, while the more nearby variation demonstrates that a 19 km s" 1 Mpc" 1 
is certainly significant, the more distant variation can not be confirmed to be 
significant. 

Unfortunately, there does not appear to be a reasonable way to use the 1- 
cr errors from the combined data sets to properly determine the statistical 
significance as in Sect. 3.2, as the errors reported for the comparison data are 
reported inconsistently and do not share the same systematic errors. However, 
just using the 36 type la supernova values, since they have small random errors 
and will share systematic errors, and since they tend to be the most distant 
values, yields interesting results. It is found that as much variation occurs 
only 0.00% of the time for a standard deviation of 35° (which yields a range 
of variation of 6.7 km s _1 Mpc -1 ), 0.02% of the time for a standard deviation 
of 25° (which yields a range of variation of 9.7 km s" 1 Mpc -1 ), and 0.12% 
of the time for a standard deviation of 15° (which yields a range of variation 
of 12.9 km s" 1 Mpc" 1 ). If the 1-cr errors are truly accurate for the type la 
supernova measurements, it means it is essentially impossible to achieve this 
much variation by chance. To the extent that these 1-cr errors can be trusted, 
it suggests the variation seen on large scales is statistically significant. 
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Fig. 4. Hubble constant contour map of Fig. 1(a) transformed to (a) celestial co-or- 
dinates, (b) ecliptic co-ordinates, (c) supergalactic co-ordinates, and (d) the CMB 
frame. The poles of the CMB frame are defined by the dipole in the CMB in the 
heliocentric frame, and the longitude right to left is defined such that the North 
Celestial Pole is at 90° longitude. Positions of the actual data points are indicated 
by triangles, and the contours range from low (dark) to high (light) values of Hq 
(in km s" 1 Mpc -1 ) as indicated. 

4 Discussion 



The variation does not appear to be an artifact of Galactic dust, since there is 
no consistent difference looking in or out of the plane of the Galaxy. In fact, the 
overall structure in the map is inc onsist ent with the distribution of dust in the 
COBE dust maps (Schlegel et al.l. ll998h . so it seems unlikely that the observed 



variation in Hq could be due to poor corrections for dust in our own Galaxy. In 
Fig. 4 the HST Key Project map of Fig. 1(a) has been transformed respectively 
to celestial co-ordinates, ecliptic co-ordinates, supergalactic co-ordinates, and 
a CMB-dipole-oriented frame of reference for comparison. Figures 4(a) and 
4(b) demonstrate that the extrema do not appear to be an artifact of any 
local frame of reference. 

The supergalactic map is interesting because the extrema that predominate 
locally are near the supergalactic plane, suggesting they may be associated 
with local conglomerations of matter (which tend to be arranged along the 
supergalactic plane). Meanwhile, the extrema that dominate farther out are 
oriented closer to the supergalactic poles, so if these extrema are a real phe- 



19 




Fig. 5. Hubble constant contour map of Fig. 2 with recent determinations of bulk 
flows ( triang les). From hig h to low latitu de, th e bulk flows are those of da Costa 
et all d2OO0h . Dekel et al| (119991 1. WilliclJ 4l999j), Hudson et al| tooi ). Parnovsky 



et al] ( 2001 ). and Giovanelli et alj (jl998l ). Positions of the minima and maxima 



are indicated by dark and light dots respectively, and the contours range from low 
(dark) to high (light) values of Hq (in km s _1 Mpc -1 ) as indicated. 

nomenon, we may be getting a clear view looking out of the supergalactic 
plane and observing effects associated with a much larger scale. 

The CMB-frame map is oriented with its north pole in the direction of the 
Sun's peculiar velocity with respect to the CMB. From the map it is apparent 
that there is little difference between directions oriented along the CMB dipole. 
All of the extrema are near the equator of the map, roughly perpendicular to 
the CMB dipole. Since the Hq values have been corrected to be in the CMB 
frame of reference, this suggests there is no error due to the correction for the 
CMB dipole: it seems too good in fact. Interestingly, this variation se ems t o 
be in agreement with the CMB anisotropy observed by Tegmark et alj (12003 ). 



who found that the CMB quadrupoie and octupole are aligned such that the 
extrema are in a plane roughly perpendicular to the direction of the dipole. 
If our dipole motion with respect to the CMB is related to the existence 
of the if /CMB anisotropy perpendicular to this direction, it suggests this 
motion and the variations related to it m ay both be stemming from large-scale 
inhomogeneities. Inoue and SillJ (j2006 ) have suggested the CMB quadrupoie 



and octupole alignment can be explained by a pair of voids a few hundred 
Mpc distant in the direction (I = 330°, b = —30°), which is interestingly in 
the direction of the H maxima in our maps. 

In Fig. 5 recent bulk flow directions are plotted on top of the map of Fig. 2. 
The bulk flow directions tend to lie on the outskirts of the uncertainty in the 
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maxima, so they are not totally consistent with the map. They should be, as 
bulk flows will be oriented toward directions associated with higher recessional 
velocities in the CMB frame and hence the directions of higher H . However, 
the bulk flow directions do lie in the higher H regions of the map, so they 
are not that disconcordant. One thing to note is that the bulk flow directions 
appear to be sandwiched between the primary and secondary maxima. This 
suggests that by only looking at the net flow, bulk flow studies may be missing 
the distinction between two separate effects and missing the actual directions 
of interest. 



5 Conclusion 

It appears that a statistically significant variation in H of at least 9 km s -1 
Mpc -1 exists in the HST Key Project data. The approximate directional un- 
certainty is 10° to 20°. Maps weighted for distance appear to indicate two 
sets of extrema that dominate on different distance scales. Within our super- 
cluster, differences as great as ~35 km s" 1 Mpc -1 are observed, and these 
tend to occur near the supergalactic plane with a minimum near (a = 9 h 30 m , 
5 = +70°) and a maximum near (a = l^SO" 1 , S = —70°). Beyond our super- 
cluster, differences as great as ~20 km s -1 Mpc -1 are observed, and these tend 
to occur away from the supergalactic plane with a minimum near (a = 18 h m , 
S = +15°) and a maximum near (a = 5 h 30 m , S = +5°). Within 70 Mpc, a 
combination of the HST Key Project data and the comparison data shows a 
statistically significant difference of 19 km s -1 Mpc -1 . Beyond 50 Mpc, the 
HST Key Project type la supernova data yield a statistically significant dif- 
ference of 13 km s _1 Mpc -1 (assuming the reported 1-cr errors are reliable). 

Further study of the resilience of this result requires more data specifically 
selected for achieving optimal sky coverage. It would also be interesting to 
have data for a greater range of distances to see how far out a statistically 
significant variation can be detected and over how large a scale the Universe's 
expansion needs to be sampled before it appears to become uniform. 

Real variation in H is not really unexpected given the degree of structure and 
mass inhomogeneity present in the Universe. One implication of this variation 
is it also partially explains why H has historically been plagued by so much 
uncertainty. 
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Abraham, Byron Desnoyers Winmill, John Dubinski, Shoko Sakai, John Tonry, 
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